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Abstract 

The first order phase transition of the three-dimensional Blume Capel 
are investigated using cooling algorithm which improved from Creutz Cellu- 
lar Automaton for the D/ J = 2.9 parameter value in the first order phase 
transition region. The analysis of the data using the finite-size effect and the 
histogram technique indicate that the magnetic susceptibility maxima and 
the specific heat maxima increase with the system volume [L'^) aX D / J = 2.9. 
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1 . Introduction 

The Hamiltonian of the Blume-Capel model is given by, 

Hj = -jJ2 SiSj + Dj2Si (1) 

<ij> i 

where Sj = —1,0,1 and the first sum is carried out over all nearest-neighboring 
(nn) spin pairs on a three-dimensional simple cubic lattice. The parameters 
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of the J and D are the bihnear interaction energy and single-ion anisotropy 
constant, respectively. 

The model is known to have a rich critical behavior and considered by 
both numerical and analytical methods. Recently, most of these studies 
mainly focused on determining the tricritical point of the model . In pre- 
vious paper*^^\ we obtained the tricritical point value of the 3-d Blume-Capel 
model as D/J = 2.82. This results agrees with a series expansion ^^~'^\ the 
cluster variation method the effective field theory and Bethe- Peierls 
approximation results^^^ . Although the critical behavior on the second-order 
phase transition region {D/J < 2.82) has been studied commonly, the finite- 
size effects in the first-order phase transition region {D/J > 2.82) is not 
investigated exactly. 

Our aim of this paper is to identify the phase transition in the first order 
phase transition region of the Blume Capel model by the finite-size effects and 
the histograms of energy distribution P(E) and order parameter distribution 
P(M). 

It is well known that the first order phase transitions involve the coex- 
istence of two distinct phases and are characterized by the existence of a 
discontinuity in the energy and magnetization for the infinite systems. As a 
results of these discontinuities, the specific heat and the susceptibility show 
singularities of the delta function. However, the characteristic singularities 
and discontinuities in the first order phase transitions appears as rounded 
and smeared in finite systems employed in computer simulations*^^^'^^). These 
behaviors at the first order phase transitions in finite-size systems are qualita- 
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tively similar to the finite-size effects at second-order phase transitions. The 
finite size effects at the first order phase transition have been investigated by 
the different methods such as the Monte Carlo'-^^"^^-*, the phenomenological 
renormahzation group^^^^ , transfer matrix method*^^^-* , renormahzation-group 
analysis^^^) . According to the results of these studies, the finite-size effects 
at the first order phase transition depend on the volume of the system L"^. 
The specific heat and the susceptibility increase with L'^ and the transition 
temperature which are locations of their extrema, approach the transition 
temperature in the infinite system as L"*^. In addition, the histograms of 
energy distribution P(E) and the order parameter distribution P(M) are 
the reliable method to study first-order phase transitions*^^^'^^'^^'^^''. While 
the probability distribution of energy shows a single peak in a second-order 
phase transition, it shows a double peak in a first-order phase transition. In 
this paper we simulate the three dimensional Blume-Capel model with the 
cooling algorithm which improve from the Creutz Cellular Automaton in the 
first order phase transition region. The finite-size scaling and the probability 
distribution for energy and order parameter are used to determine the nature 
of the phase transition. The remainder of the paper is organized as follows. 
The data are analyzed and the results are discussed in Section 2 and the 
conclusion is given in Section 3. 
3. Results and discussion 

The three dimensional Blume-Capel model is simulated with the cool- 
ing algorithm*^^^ which improved from the Creutz Cellular Automaton. The 
simulations have been made at the D/J — 2.9 anisotropy parameter value 
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in the first order phase transition region on simple cubic lattices LxLxL of 
the linear dimensions L — 8, 10, 12, 14, 16, 18 and 20 with periodic boundary 
conditions. 

At the algorithm, the cooling rate is equal to O.OlHk per site for D/ J — 
29/10 value and the kinetic energy of the system reduced by the different 
cooling amounts per site because the kinetic energy, if^, is an integer variable 
in the interval (0, 24 J) . The computed values of the quantities are averages 
over the lattice and over the number of time steps (1.000.000) with discard 
of the first 100.000 time steps during which the cellular automaton develops. 

The simulations were done about 10 times with different initial config- 
urations at each lattice size. Measurements shown that the variation of 
thermodynamic quantities and the histogram of P(E) and P(M) for each 
different initial configuration are different from each other around the transi- 
tion temperature. At the same time, some of the simulations do not exhibit 
the double peak in the histogram of the P(E) and the three peaks in the 
histogram of the P(M) which characterize the first order phase transition. 
Therefore, we choose the simulations which exhibit the characteristic double 
peak in the histogram of P(E) and the three peaks in the histogram of the 
P(M) to study finite-size effects at D/J = 2.9. 

In Fig.l, the temperature variation of the order parameter (M = -p- X]<ij> 
and Q = E<y> Sf), the energy {E ^ -JJ2<ij> SiSj/E^ where Eq is the 
ground state Ising energy at kT/J=0), and histograms of P[M) and P{E) 
are shown for a chosen simulation result on L = 18 as an example. There are 
discontinuities which characterize a first order phase transition in tempera- 
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ture dependence of the order parameter and energy in Fig.l (a) and (c). On 
the other hand, another evidence for the first order transition is seen in the 
histograms of P(M) and P(E) in Fig.l (b) and (d). The shape of distribution 
P(E) shows a double peak around the transition temperature corresponding 
to the coexistence of the ordered and disordered phases. In the distribution 
of the magnetization P(M), the middle peak indicate to the value of disor- 
dered phase while the other peaks indicate the order phase. The histograms 
of P(M) and P(E) around the transition temperature is given in Fig.2(a) and 
(b) for a chosen simulation, respectively. Histograms of P(M) and P(E) have 
a single peak above the transition temperature. As the temperature is close 
to the transition value the histogram of P(E) has a double peak while P(M) 
has three peaks indicate the coexistence region. With the decreasing tem- 
perature, the histogram of P(E) converges to a single peak and the middle 
peak in the histogram of P(M) is depressed while other peaks are enhanced 
corresponding to order region. The behavior of P(M) and P(E) around the 
transition region are in good agreement with the characteristic behavior of 
the first order phase transition (I0'i2,i3,i9) 

In Fig. 3 (a) and (b) the temperature dependence of the order parameter, 
the magnetic susceptibility, the energy and the specific heat are shown at 
different lattice sizes for a selected simulation which shows a double peak in 
the histogram of P(E) and three peaks in the P(M). The results are seen 
in Fig. 3 show that the energy and the order parameter have discontinuity 
at transition region and the peaks of the susceptibility and the specific heat 
increase with lattice size as expected in the first order phase transition. 
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The finite-size effects of tlie specific fieat and the susceptibihty at the 
transition temperature are given by 

C'max Of L'^ (2) 
Xma. « L'^ (3) 

The magnetic susceptibility maxima Xmaxi.^)^ the specific heat maxima 
Cmax{L) are obtained from the average over the maximum of magnetic sus- 
ceptibility and specific heat at each lattice size for chosen simulations. How- 
ever, the transition temperature Tt{L) for each lattice is estimated from the 
location of the magnetic susceptibility and the specific heat maxima for cho- 
sen simulations. The infinite lattice transition temperature are estimated 
from the extrapolation of susceptibility T^{L) and specific heat maxima 
Tf(L) for various lattice sizes. The estimated infinite lattice transition tem- 
perature vahics arc T^{oo) = 1.29±0.02, Tf{oo) = 1.28±0.03 from Xmax 
Cmax, respectively. The logarithm of the specific heat and the susceptibility 
maxima as a function of the logarithm L is shown in Fig.4(a) and (b). The 
data lies on a single curve and the slope gives 3.01 and 3.07 from the specific 
heat maxima Cmax{L) and the susceptibility maxima Xmaxi^), respectively 
These estimated values are in good agreement with the system dimension 
{d — 3). This result show that, the magnetic susceptibihty and the spe- 
cific heat maxima increase with the system volume (L*^) as a characteristic 
property of the first order phase transition 

4. Conclusion 
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The three dimensional Blume-Capel model is simulated using cooling al- 
gorithms on a cellular automaton. The simulations are done about 10 times 
with different initial configurations at each lattice size and some of this sim- 
ulations having double and three peaks in the histograms of P(E) and P(M) 
are chosen. To identify the phase transition of the model, the analysis of 
the first order phase transition effects and the histogram technique are used. 
The temperature variations of the order parameter and the energy show the 
discontinuity at the transition temperature while the histogram of P(E) and 
P(M) have a double and three peak structures. Furthermore, the magnetic 
susceptibility and the specific heat maxima increase with the system volume 
{L'^) as expected in the first order phase transition. 
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Figure Captions 

Fig.l. The temperature dependence of (a) the order parameters M and 
Q, (b) the histogram of P(M) ,(c) the energy, (d) the histogram of P(E) for 
a chosen simulation on L=18. 

Fig. 2. The histogram (a) P(M) and (b) P(E) for several temperatures 
around the transition temperature for choosen simulations on L =10 and L 
= 18. 

Fig. 3. The temperature dependence of (a) the order parameters, (b) 
magnetic susceptibility, (c) the energy and (d) the specific heat on L = 
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8,10,12,14,16,18 and 20. 

Fig.4. Log-log plot of (a) the specific heat maxima and (b) the magnetic 
susceptibility maxima against linear dimension L. 
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